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FLOW THROUGH HORIZONTAL VENTS 
AS RELATED TO COMPARTMENT FIRE 
ENVIRONMENTS 


QING TAN 
AND 
YOGESH JALURIA 


ABSTRACT 


A detailed investigation has been carried out on the flow exchange through a horizontal 
vent in a compartment containing a fire. A plexiglass tank with a vented horizontal 
partition in the middle was constructed to simulate the warmer interior environment 
due to a fire and the cooler ambient environment by filling the upper and lower com- 
partment with brine and pure water, respectively. Experiments have been carried out 
on the combined natural and forced convection flow by imposing a pressure difference 
across the vent. The flow rates through the vent were determined over wide ranges of 
the governing variables, such as the pressure difference AP across the opening, density 
difference Ap across the opening and the opening length to diameter ratio L/D. The ba- 
sic characteristics of the flow, particularly whether it is unidirectional or bidirectional, 
was also studied. Volume flow rates were obtained as functions of the governing pa- 
rameters in terms of correlating equations, from which quantitative information of the 
effect of AP, Ap and L/D on the flow exchange through the vent can be determined. 
These results can thus be applied to the modeling of fire growth in vented rooms. 
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Chapter 1 


Introduction 


1.1 The Problem 


The study of the flow of the combustion products out of a room with a fire and the mass 
flow rate of oxygen being added to the room by the flows at the openings and the vents 
of the room is of considerable importance in the mathematical modeling of the fire and 
of the changing environment in the room[24]. It also has a number of applications in 
many other analogous situations. One application involves the flow that arises during 
the venting procedure employed by firefighters conte: combating room or building fires. 
Another application concerns fluid motion through apertures connecting various rooms 
in buildings. Street drains, partially filled with water that is warmer than the ambient 
air, represent yet another application. This is especially important if a highly volatile 
toxic liquid or gas that is lighter than air finds its way into the drainage system [17]. 
Similarly, many environmental and geothermal flow processes, flows in nuclear reactor 
vessel systems and flows in passive solar heating energy systems give rise to vent flows 


that are often of crucial importance in such systems. 


Among various vent flows, the exchange flow through a small vent opening in a 
horizontal partition is of particularly importance and needs to be studied in detail. 
The rate of movement of radiological or toxic gases in tall industrial or commercial 
buildings under emergency conditions is mainly a function of the leakage between the 
floors rather than that through the walls. Also, from a fundamental point of view, there 
is a big difference between horizontal and vertical opening or passage with respect to 
the relative motion of the two fluids. The flow pattern in a vertical passage through 


a horizontal partition, such as ceiling or floor vent, is often disturbed and chaotic and 


may not be amenable to analysis. However, the two fluid regions in a horizontal pas- 
sage through a vertical partition have been observed to be stably stratified and easily 
analyzed by Bernoulli’s equation [19,22]. 


Basically, the flow exchange through a horizontal vent involves a combination of a 
buoyancy-driven exchange flow and a forced flow . The density difference Ap across 
the opening between the cold (heavier) ambient fluid and the warm (lighter) enclosure 
fluid arises partly from the difference in composition and partly from the difference in 
temperature due to the heating up of the gases in the enclosure containing the fire. The 
density difference gives rise to a buoyancy-driven downflow of the heavier fluid from the 
ambient environment into the enclosure. However, if the compartment is sealed and has 
no other openings, mass conservation dictates an equal volume upflow of the lighter gas. | 
The pressure difference AP across the opening arises from a momentum force balance 
in the flow, with the pressure in the enclosure higher than outside. The pressure differ- 
ence generates a forced upflow of the lighter gas against the buoyancy-driven downflow 
of the heavier ambient fluid. The resulting flow rate through the vent, as a function of 
Ap, AP and the thickness to diameter ratio L/D of the vent, is of direct relevance to 
fire research. Here AP = Pgor - Prop, and Ap = prop - psBor, the subscripts “TOP” 
and “BOT” refer to the regions above and below the vent, respectively(Fig.1.1). Thus, 


positive values of AP and Ap are of particular interest in typical enclosure fires. 


Some work has been done on buoyancy-driven vent flows for the special case of AP 
= 0.0. Work has also been done on flows due to an applied pressure difference, with no 
buoyancy effects present. But very little information is available for nonzero pressure 
differences combined with buoyancy effects, this circumstance is of particular interest 
in vented room fires. The flow exchange needs to be quantified over wider ranges of 


Ap, AP and L/D ratio in order to provide the necessary inputs for fire modeling studies. 


Figure 1.1: The basic horizontal vent configuration. 


1.2 Literature Survey 


Buoyancy-driven flows in enclosures have received a tremendous amount of atten- 
tion[{1,2]. However, not much work has been done on the flow through openings, such 
as those between connecting rooms in buildings. Vertical openings have received some 
attention because of their importance in several practical problems, such as electronic 
and energy systems and room fires. Fig.1.2 shows a sketch of the frequently consid- 
ered circumstance of fire in a compartment with a single vertical opening (a door or a 
window). In this model, the hot gases rising above the fire, toward the ceiling, from a 
distinct layer of the hot gases. Air flows into the room at a rate mg, shears or entrains 
some gas from the hot upper layer at a rate m,;, and mixes with it in the lower layer. 
The fire plume entrains gas from the lower layer and “pumps” it to the upper layer at 
a rate m,.. Under steady conditions the mass flow rate through the opening from the 


upper layer (m,)equals the fuel flow rate plus 7,[8]. 


Experiments on the exchange flow through openings (windows or doors) of small L/S 
ratio, which is the ratio of opening length to the side of square opening, were reported 
by Brown and Solvason [3] and Shaw [4]. Shaw also studied combined buoyancy-driven 
and forced air flow across a door in a vertical partition. Buoyancy-driven flow through a 
horizontal duct or tube, in a vertical partition, between two compartments was studied 
experimentally by Leach and Thompson[5] and Bejan and Rossie[6]. Experimental and 
theoretical studies of steady state, fire-induced flow through vertical openings were also 
conducted by Prahl and Emmons [7], Steckler et. al. [8], [9] and Emmons/10]. In the 
later paper by Steckler et. al. [9], the brine-water technique was used to model the 
buoyancy driven flow of combustion products within a multicompartment structure. 
Several numerical studies have been carried out on this problem, such as those by Chan 
and Tien [11] and Abib and Jaluria [12]. As mentioned before, the flow across a vertical 
opening is stably stratified and can be obtained by using Berioulli’s equation and a jet 


contraction coefficient to estimate the velocity distribution at the vent. 
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Figure 1.2: Fire in a compartment with a single vertical opening 


The work done on the flow through horizontal vents, such as the one shown in 
Fig.1.3, is very limited. From a survey of the literature, it appears that the only prior 
papers dealing with vertical buoyancy-driven exchange flow across a horizontal parti- 
tion between two compartments connected by a single opening are those of Brown[13] 
and Mercer and Thompson [14]. In Brown’s experiments, air was used as the fluid 
and the countercurrent flow through a square aperture was investigated by imposing 
a temperature difference between the two compartments. His test covered ratio L/D, 
of the opening length(or partition thickness) to the side of the square opening, in the 
range of 0.0825 to 0.66. Brown interpreted the countercurrent flow as a heat transfer 
phenomenon and expressed his results in terms of a correlation for the Nusselt number 
Nu, versus Grashof number Gr. By introducing the definitions for Nu and Gr, and 
canceling viscosity and thermal conductivity, his correlation can be re-expressed, for 
practical purposes, in the form of a single functional relationship between the exchange 
flow rate and Ap/f and L/S. Here, Ap is the density difference across the aperture and 
p is the mean density between two compartments. A major result of Brown’s investi- 


gation is that the exchange flow rate increases with increasing L/S. 


Mercer and Thompson’s [14] experiments investigated the buoyancy-driven flow 
through inclined short tubes by using brine and water as the working fluids to ac- 
complish the density difference between two compartments. Their experiments with 
vertical tubes, which are of interest here, covered a tube length-to-diameter ratio(L/D) 
range of 3.5 to 18.0. Interestingly enough, the exchange flow rates were observed to 
decrease with increasing L/D which is opposite to that observed by Brown. They sug- 
gested that, for a single opening, a maximum exchange flow rate exists at L/D or L/S 


between 0.66 and 3.5. 


Several other investigations have been carried out on the buoyancy-driven flow 
through a horizontal vent. These include papers by Thomas et al. [15], Hinkley[16], 
Myrum[17], Cooper[18], Epstein[19], and Mercer and Thompson [20]. All of these stud- 


ies mainly concentrated on the special circumstance of AP = 0.0. Air was used as the 
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Figure 1.3: Fire in a compartment with a single horizontal opening 


working fluid in some of the references while water/brine for the rest. 


Employing water/brine as the working fluid, Epstein [19] experimentally considered 
the effect of varying the L/D ratio on the exchange flow rate and identified four different 
flow regimes. Fig.1.4 gives a sketch of the experimental arrangement used and the major 
results obtained. An empirical correlating equation for the dimensionless exchange flow 


rate, in terms of a Froude number Fr, versus L/D ratio is expressed as below: 


SS 2Q) ae 0.055[1 + 400(L/D)3}*/6 
(ghpD®/ayi2 ~ [1 +0.00527(1 + 400(L/D)5)¥/2((L/D)° + 117(L/D))5/4}/5 
(1.1) 


Lr 
The parameter ranges covered are : 


0.01 < L/D < 20.0 and 0.025 < Ap /p < 0.17 


where Q is the volume exchange flow rate, g the magnitude of gravitational accelera- 


tion, Ap the density difference across the vent and p the mean density. 


It may be mentioned here that the characteristic velocity Vo in a buoyancy-induced 


flow is given by[1,2] as: 


Vo = VgApD/p (1.2) 


since the area of a circular vent is xD*/4 , the flow rate varies as V-rD?/4. Thus, this 


quantity is used to nondimensionalize Q, yielding the Froude number. 


The various flow regimes were termed by Epstein as Oscillatory Exchange Flow, 
Bernoulli Flow, Combined Turbulent Diffusion and Bernoulli Flow, and Turbulent Dif- 
fusion Flow. The results of Epstein’s experiments suggested that there exists an in- 
fluence of the density ratio on the Froude number Q/(gApD*/p)'/2, especially at low 


L/D where an increase from Ap/p = 0.04 to 0.09 was observed to give rise to modest 
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Figure 1.4: Sketch of Epstein’s experimental arrangement and his major results 
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reductions in the Froude number (~ 20 percent). A further increase in Ap/p above 0.09 
or, for that matter, decrease below 0.04, were found not to lead to significant changes 
in the Froude number. Thus, it was concluded that no additional density parameter is 
needed to characterize the results. The flow rate attains a maximum value at an L/D of 
about 0.6. At lower L/D the flow rate decrease with decreasing L/D and appears to be 
tending toward a constant value as L/D -0. At the higher L/D the opposite behavior 
is noted, with the flow rate decreasing with increasing L/D. This might explain the 


conflict between the results of Brown [13] and Mercer and Thompson [14]. 


The conditions of steady, purging flow between space connected by a relatively 
long, straight, inclined, circular pipe of diameter D have been determined experi- 
mentally by Mercer and Thompson [20]. The configuration for their experiments is 
sketched in Fig.1.5, The bottom and top tanks contained fresh water and a relatively- 
dense salt-water solution, respectively . Each experimental run involved a steady 
uni-directional flow of salt-water, from the outlet of the top tank, at a measured 
volume-flow-rate, Qayeas. For small enough Qyszas and at the entrance to the 
lower tank, a wedge of the lower-tank fresh water was observed flowing back into the 
pipe. The length of the flow wedge W increased with decreasing Qyezas and W/D 
was taken as a measure of the purging effectiveness associated with the dimensionless 
flow rate, Qarras/(gApD*/prop)/*. A total purging implies a unidirectional flow 
through the vent from one space to the other, depending on the sign of AP, where 
AP = Pgor— Prop. A minimum totally purging value of Qyzas is defined here as a 
flow rate which would lead to W = 0 exactly. The “purging” terminology in the above 
discussion is taken from the technology for nuclear reactor safety, where it is important 
to be able to predict conditions where totally purging flows will give way to flows involv- 
ing nonzero exchange flows between reactor vessels and outside fluid sources. As seen 
in Fig.1.6, the W/D vs. Qucas/(gApD*/prop)'/? data are a strong function of pipe 
inclination, 6, where 8 = 90 (vertical pipe in a horizontal partition) is of special interest 
in the present study. It is impossible to observe the wedge length for values of W/D 


less than about two because of oscillations. The critical purging flow rate is therefore 


yak 


Figure 1.5: Configuration of Mercer and Thompson’s experimental set up 
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Figure 1.6: Experimental results of Mercer and Thompson’s: Non-dimensional wedge 
length as a function of the non-dimensional velocity: results for 0°, 45°, 60°, and 90° 
inclinations. 
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taken as that corresponding to a wedge length-to-diameter ratio of two. Fig.1.7 shows 
the purging flow rates required to give a minor W/D ratio (W/D = 2) as a function of 
Ap/ p and inclination 6. 


Few studies have considered the combined natural and forced convection through 
horizontal vents. Cooper[21] calculated the flow through a horizontal vent located in a 
ceiling or floor of a multi-room compartment. he took the cross-vent pressure difference 
AP into consideration in addition to the density difference. He assumed that in each 
space (bottom or top), near the vent elevation, but away from the vent flow itself, the 
environment is relatively quiescent with the pressure well-approximated by the uniform 
pressure field which would be derived from hydrostatic momentum considerations, that 
is tosay AP = Pgor-Prop, where Pgor and Prop are the hydrostatic pressures at the 
vent elevation. By employing the experimental results of Epstein [19] and Mercer and 
Thompson [20] and the Standard Vent Flow Model[10], Cooper developed an algorithm 


for flow through horizontal vents. 

The Standard Vent Flow Model, based on Bernoullis’s equation, assumes unidirec- 
tional flow and breaks down if a bidirectional flow exchange occurs under the effects of 
buoyancy and the applied pressure difference. According to the Standard Vent Flow 
Model, for different AP and with prop > pot, the volume flow rate Q s07,sT/QrTopPsT 
through the vent, as illustrated in Fig.1.8, is given by: 


For arbitrary Ap = prop — pBoT: 


Qsorst = CpAy|24P/pgpor|'’?, Qropsr=0 ifAP>0 (1.3) 


Qso7,sT = 0, Qropst = 0, if AP = 0 (1.4) 


Qsorst=0, Qropst=CpAv|2AP/prop|?, if AP <0 (1.5) 


1 /min 


Purging flow ; 
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Figure 1.7: Experimental results of Mercer and Thompson’s: Full-scale purging flow 
requirements for a wedge length of twice the diameter(D = 1.8 m) 
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(a) AP>0: Q Topst =0 


(b) A P= 0:Q Top,st= Q Borst =0 


(c) 4 P<O : Q BOT,ST=0 


Figure 1.8: The standard vent-flow model for horizontal vents. 
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where Qgor,sT and Qropsr are the volume-flow-rates through the vent from the bot- 
tom space to the top space and from the top space to the bottom space, respectively, 
subscript “ST” stands for the Standard Vent Flow Model, Ay is the vent area, and 
C'p is the flow coefficient. So, the Standard Vent Flow Model predicts a zero flow rate 
when AP = 0, which is obviously wrong when a heavier fluid overlies a lighter fluid, 
because of the resulting instability, it is evident that a two-directional or an exchange 


flow, develops where the top and bottom spaces exchange fluid across the opening. 


By assuming the correct vent flow, Qrop and Q gor, can be predicted by a sum of 
the standard model flow components of Eqs.{3] - [5] and an exchange-flow component, 


Qzx, Which is to be determined. Cooper[21] assumed that: 


Qtop = Qropsr + QEx (1.6) 


Q sot = QsBor,st + QEx (1.7) 


where he further assumed that Q x is a linear function of AP/|AP oop], as shown in- 


Fig.1.9. Then from Epstein’s experimental results without pressure difference: 
Qex = Qex,max = 0.055[gApD*/p]"/”_ at AP = 0 (1.8) 


The pressure at which Qzx = 0 is the flooding pressure APrroop, which is the 
critical point leading to the totally purging flow from one side of the vent to the other, 
depending on the sign of AP. Again, by assuming that the exchange flow through the 
vent is irrelevant to vent inclination when it is ae close to the totally purging flow, 
Cooper deduced the value of APp,o0p based on Mercer and Thompson’s experimental 


results as below: 


|APrroop| = CraaprgApD*/ (2A, (1.9) 


where Csyapp is the shape function. So, when |AP| = 0, there is the maximum coun- 


tercurrent exchange flow due to the buoyancy force. As |AP| increases from zero, but 
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Q/(294pD*/(erop + pBor))*/? 


Q=Qrop ST 


AP/|APrtoop!| 


Figure 1.9: 'The volume-flow rate components Qropsr and Qzx for unstable Tig.1.8 
configuration, Ap > 0, plotted as a function of AP. 
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is less than |APpyoop|, there is some flow exchanges, i.e., non-zero Vex, between the 
environments of the two spaces. If |AP| keeps increasing and exceeds |A Pproop|, there 
is no exchange flow and the flow through the vent is uni-directional. Following the 
Standard Vent Flow Model, this leads to a total purging flow from one space to the 
other determined by the sign of AP. The VENTCL algorithm developed by Cooper is 


suitable for general use in zone-type compartment fire models. 


Epstein and Kenton[22] investigated the combined natural convention and forced 
flow through small openings in a horizontal partition based on the former research of 
Epstein{19]. Two different experimental apparatus had been employed by them. One 
was with the draining technique and the other with the water injection technique. These 
are shown in Fig.1.10 and Fig.1.11, respectively. Separate experiments, with these two 
methods, were performed on forced flow through the vent opening to obtain a definitive 
relation between the flooding flow rate q, above which the buoyancy-driven exchange 
flow is prevented, and the opening length-to-diameter ratio L/D. An empirical relation 
for the flooding (purging) flow rate in terms of the Froude number versus L/D was 
found as depicted in Fig.1.12 and the following correlation was given with a relative 


error of less than about 20 percent: 


q 0.19[1 + 4 x 109(L/D)5}?/9 


(gApD®/py? ~ (145.091 X10-°(L/DP 71 + 4x 10L/ Dey (1-10) 


Here q is the purging volume flow rate. Interestingly enough, the shape of the curve 
is similar to the exchange flow rate versus L/D trend reported in Epstein’s previous 
paper, which was without the presence of an external pressure. The main intention of 
their later study was to determine the effects of forced-fluid motion on natural convec- 
tion for horizontal venting flow. A correlation for this flow was obtained and is shown 
in Fig.1.13. Here Qgr is the buoyancy-driven component of the combined convection, 
Qcc the purely buoyancy-driven exchange flow rate with AP = 0.0, Qy the net sup- 


plying flow rate and q the purging flow rate. It can be seen from the Fig.13 that the 
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Figure 1.10: Schematic diagram of the Epstein and Kenton’s apparatus for the deter- 
mination of the flooding flow rate via a quasi-steady draining technique. 


20 


worz2efut 1a3eM Aq 23838 Apeazs 48 94e1 MOQ ZUTpPooy ay} Jo UOTZeUTUI 
~19}2P 2q} Joy snjeredde s.uojuay pue urezsdq ay} Jo UIVIZEIP 91QeUIAYIG :[ TT eM3ry 


“ES __HaLVM 

2: ee =e HS3uY4 
YALUM A 

b ANIVA 

ef aeal face 

i \ 

(kw! - 
>ANV 1 INIUG 


MO 154SAO 


y) 
\ 


D,S (m) 

© 0.0254 

CO) 0.0445 

BM 0.0292 (Square) 


q 
(D°gApy/p)'/? 


0.01 0.1 1.0 10.0 


£ 
D 


Figure 1.12: Epstein and Kenton’s experimental results for the critical purging (or 
flooding) flow rate into a single opening in a horizontal partition 
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Figure 1.13: Epstein and Kenton’s experimental results: The natural convection com- 
ponent Qgr of combined natural convection and forced flow across circular openings 
in a horizontal partition; correlation of data in terms of Qcc, Qu, and q. 
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externally supplied flow rate Qy need only exceed a value equivalent to about 60 per- 
cent of the flooding rate q to essentially “wipe away” the natural convection component 
Qpr. An expression for the solid, correlating line in the Fig.13, which gives Qgpr to an 


accuracy of 25 percent, is: 


Qar = Qcc(1- Qu/q)* (1.11) 


this equation compares well with the theoretical prediction. However, only the volume 


flow rate value was determined and no pressure information was given. 


1.3 Objectives 


As reviewed in the preceding sections, the flow through horizontal vents, such as those 
in the ceiling or the floor of a multi-room compartment, is important in the modeling 
of fire in a vented enclosure. Such flows are driven by buoyancy effects, due to the 
difference in density of the fluids on either side of the vent , as well as by any pressure 
difference that exists across the vent. Some work has been done on the flow through 
horizontal vents due to the density difference for a zero pressure difference. Similarly, 
the Standard Vent Flow Model works well for unidirectional flow and for circumstances 
where the density difference is small or the distribution is stable. Some experimental 
work has also been done to study the combined effect of AP and Ap and to quantify 
the purging effectiveness of the flow. Clearly, the flow through horizontal vents is an 
important consideration in the modeling of compartment fires for several circumstances 
of practical interest. Due to the lack of detailed, quantitative information on the flows 
exchange as a function of arbitrary values of the governing variables AP, Ap and L/D, 
the main objective of this work is to carry out a detailed investigation on this combined 
natural and forced convection phenomena. The flow rate through the horizontal vent 
was determined over wide range of the governing variables. The basic characteristics 
of the flow, particularly whether it is unidirectional or not, was also studied. The main 


thrust of the study was to obtain the flow exchange rate through a horizontal vent as a 
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function of AP, Ap and L/D and to present the results obtained in terms of correlating 


equations which may by readily applied to modeling of fire in vented rooms. 


The experiment made use of the fact that saline water is denser than fresh water. 
This attractively simple choice of fluids to represent buoyancy-driven gas movement 
between rooms is valid so long as molecular viscosity and thermal diffusivity are not 
important parameters. The brine-water technique has some advantages over performing 
a heat transfer experiment with gas as the working fluid. There is no need to cover the 
test section with insulation as the flow exchange process is isothermal, and there is no 
significant waiting time requirement for the steady state to be reached. The maximum 
compartment-to-compartment density difference ratio Ap/p that can be achieved with 
the brine-water system is approximately 0.2. This corresponds to a gas exchange rate 
driven by a temperature difference of 100K between two rooms when the average tem- 
perature is 500K. The nondimensionalization will allow the use of the results obtained 
for modeling of air and gases, which are obviously of greater interest in fires. Since the 
water is incompressible, the application of the pressure difference AP in water is also 
relatively easy. Comparison of the small-scale Berepeentat results with previously con- 
ducted fire tests in a full-scale facility has demonstrated the utility of the brine-water 


technique in predicting full-scale results with hot gases. 


For the buoyancy-driven flow across a horizontal vent of small height L, such as the 
one shown in Fig.1.3, the velocity level may be estimated to be of order (gApD/j)'/? by 
using Bernoullis’ equation. Therefore, for a circular opening, the flow rate Q varies as 
(gApD*/p)'/?, which may be taken as the characteristic volume flow rate for nondimen- 
sionalizing Q, as done by Epstein{19]. For a vent of significant height L, the parameter 
L/D arises as an additional variable. Considering now the forced flow due to a pressure 
difference AP, Bernoulli’s equation gives an estimate for the velocity across the vent 
(2AP/p)!/?. Therefore, a buoyancy parameter B and the dimensionless flow rate Fr as 
Froude number may be defined as: 
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a (gApD®/p)1/? (1212) 


As B— 0, the forced flow circumstance, with negligible buoyancy effects, is obtained, 
and as B- oo, the pure buoyancy-driven flow, with no externally imposed pressure 
difference AP, is achieved. Thus, the flow rate Q is to be determined at different 
values of B, particularly for the mixed convective circumstance where both transport 
mechanisms are significant. 

In summary, the experimental study quantitatively determined the volume flow rate 
through horizontal vents over a range of the governing physical parameters AP, Ap and 
L/D. Employing the fresh-water /salt-water system, the flow exchange was studied both 
by visualization and by measurement of the transport across the vent. The flow rate Q | 
was obtained as a function of these variables and the results were presented as correlat- 
ing equations which can be easily applied to modeling of flows in rooms with horizontal 


vents. 
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Chapter 2 


Experimental System 


2.1 Experimental Arrangement 


Fig.2.1 and Fig.2.2 show a sketch and a photograph of the experimental arrangement. 
A rectangular tank has been fabricated with the interior region 0.373m long, 0.432m. 
wide and 0.622m high. The tank is made of plexiglass so as to allow flow visualization, 
for which a shadowgraph is used. A horizontal plexiglass partition plate located 0.311m 
above the bottom of the tank divides the tank into an upper compartment and a lower 
compartment. The lower compartment is filled with pure water and the upper one 
with brine. The partition plate is bolted with a gasket onto a shelf constructed along 
the wall. The gasket can effectively prevent leaks between two compartments. At the 
center of the partition there is an aperture whose cross-sectional area is large compared 
with the vent openings used in the experiments. A steel frame, which is made of angle 
iron, is constructed around and underneath the plexiglass tank to support and prevent 


bulging of the tank. 


Two types of openings were employed in the experiments. In one case, the brine- 
water flow was directed through a simple orifice, which was constructed by cutting a 
circular hole in a plexiglass or metal plate. In the other, the brine-water flow passed 
through a plexiglass tube that was bounded by a plexiglass support plate. The plexiglass 
orifice plate, or the tube support plate, was mounted on the partition. As illustrated 
in Fig.2.1, the centerline of the orifice or the tube was See with the centerline of 
the aperture cut in the partition. The support plate and the partition were pressed 
together by six bolts and a rubber gasket. Again, the rubber gasket was to prevent 


leaks between the compartments. Twelve round openings of different diameters and 
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(b) Back view. 


Figure 2.2: Views of the experimental set up 
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lengths were fabricated from tubes or holes cut in the support plate. The diameters of 
the openings varied from 0.009525 to 0.03175 m, and the length of the openings varied 
from 0.00127 to 0.076 m. This combination of orifices and tubes resulted in the L/D 
ratio varying from 0.05 up to 6.008. L/D < 0.25 were orifices and L for this configura- 
tion actually corresponds to the thickness of the orifices plate. The L/D > 0.25 cases 
were achieved with tubes or ducts, and L for this geometry is simply the tube or duct 
length. Tests were carried out with average brine densities over the duration of a run 
in the range of 1025-1180 kg/m°. The corresponding average density difference ratio 
range was Ap/p = 0.024 to 0.17. 


A notch is located on the left hand side wall of the tank with its base line 0.572 m 
above the bottom of the tank. This was designed to maintain the salt water level in 
the upper compartment at a constant height. The upper compartment is opened to the 
atmosphere and, therefore, the pressure is simply the hydrostatic distribution in the 
salt water. A small plexiglass container with a duct at bottom, which is designed to 
collect the overflow salt water, is attached to the notch (Fig.2.1). It works like a guide 
tubing for the notch to prevent disturbances to the overflow, and it can be closed at 
the bottom with a rubber stopper during the mixing of the salt water. Thus, it serves 


as an overflow collecting box. 


The pressure difference AP between two compartments was imposed by a translu- 
cent plastic storage tank of fresh water. The storage tank was set on a movable platform 
with a lifting jack underneath it. The bottom of the storage tank was connected to 
the bottom of the lower compartment’s right hand side wall via a tygon tubing with 
a control valve at the lower compartment end(Fig.2.1). The storage tank also has a 
notch on the side wall, so it can maintain a constant water level when fresh water enters 
the tank from a faucet. By vertically moving the storage tank up and down with the 
lifting jack to a fixed position, and with the control valve being opened, the lower com- 
partment was then pressurized according to the height of the water level in the storage 


tank. The pressure in the upper compartment is the atmospheric pressure plus the 
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hydrostatic distribution in the salt water, and the pressure in the lower compartment 
is the atmospheric pressure plus the water head of the storage tank. Thus, different 
values of AP between two compartments could be obtained. In the experiments, the 
pressure in the lower compartment was generally kept higher than that in the upper 


compartment. 


Five manometer taps are located on the same side wall of the plexiglass tank with 
the notch opening. They are all on the vertical center line of the side wall and are 
0.0127, 0.1524, 0.2731, 0.3334 and 0.4731 m above the bottom of the tank, respectively 
(Fig.2.1). Three of them are in the lower compartment while two are in the upper one. 
These are used to measure the tank pressure at different height locations, and, therefore, 
the pressure difference AP between two compartments can be deduced. By connecting 
the tap to the manometer through the tygon tubing, the pressure is obtained. The pres- 
sure readings were obtained directly from the manometer and also from the differential 
pressure transducer (Omegabrand, Model PX154-001DI). The pressure transducer was 
placed at the top of the steel frame and was on the side the closest to the manometer 
taps. The output of the transducer was connected to a multimeter (Hewlett Packard 
Brand, Model 3466A) and a chartrecorder (Omegabrand, Model 585) which were placed 


at lower front corner of the plexiglass tank. 


The shadowgraph technique was applied here for visualization purposes. The ar- 
rangement used in the shadow method, which is usually attributed to Dvorak(1880), 
is the simplest of all optical visualizing procedures[25]. As illustrated in Fig.2.3 and 
Fig.2.4, a projector was placed at the back of the tank while a smoke glass screen was 
placed in the front. The almost parallel light beam generated by the projector passes 
through the section around the vent opening. Because of the different deflection of 
light for different density fluid, a shadow of the density variation is projected onto the 
screen. One could readily observe the phenomenon or take a photograph of the shadow 
on the screen. The sharpness of the image obtained on this screen depends on the size 


of the light source. By tracing the most-separated rays emerging from an extended light 
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source, one finds that the lack of image sharpness is approximately given by Id/ f,[25], 
where d is the source diameter, / the distance between the smoke glass screen and the 
center of the test section, and f; the focal length of the lens in the projector. It follows 
that the light source should be small-not smaller, however, than a certain limit where 
the lack of sharpness due to diffraction effects becomes appreciable. The projector was 
placed 1.0m away from the tank and the screen 0.05m in order to get the clearest shad- 


ows on the screen. 


In order to mix the fluid in the lower compartment after running the experiment, a 
magnetic stirring mixer(Fisherbrand model 11-500-12SH) was employed. Unlike the up- 
per compartment, the lower compartment of the tank is hard to reach when the opening 
is closed. Lying at the bottom of the tank, the stirring bar (Magnet) was driven by the 
magnetic field generated by the mixer underneath the plexiglass tank and fulfilled the 
task very well without drilling a hole in the lower compartment(Fig.2.1). The mixing 
is much easier in the upper compartment, since it can be done manually by stirring the 


brine with a handy mixer. 


A drainage hole is located at the bottom corner of the right hand side wall of the 
plexiglass tank and is connected to a hose via a control valve. This was made to dis- 


charge the contents after the experiment. 


2.2 Experimental Procedure 


2.2.1 General Description 


To carry out the experiment, the lower compartment and the storage tank were filled 
with pure water. The opening(orifice or tube) across the partition was then closed with 
a rubber stopper and brine was introduced into the upper compartment. The pure 
water was filled up to the overflow notch of the storage tank and the brine up to that 


in the upper compartment. The storage tank was brought to a certain height with 
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Figure 2.4 
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its water head being maintained at a desired value. The control valve connecting the 
storage tank and the lower compartment was kept open, so was the opening of the duct 


in the overflow collecting box. 


A run was initiated by removing the rubber stopper. A countercurrent exchange 
flow or a unidirectional flow of brine or water was observed close to the opening by 
means of the shadowgraph. The pressures in the upper and lower compartments were 
measured. The opening was then closed, the brine density in the upper region and the 
volume overflow rate were measured. All of these were accomplished over regular time 


intervals. Further details on these measurements are given below. 


2.2.2 Density measurement 


Hydrometers(Fisher brand Precision Specific Gravity Hydrometers) have been em- 
ployed in the experiment to determine the brine density. Hydrometers of five different 
ranges, with the lowest one being from 1.000kg/dm® to 1.070kg/dm* and the highest 
one from 1.240kg/dm* to 1.310kg/dm*, were used to ensure high accuracy in the den- 
sity reading. The resolution of the lowest range hydrometer is as small as 0.1kg/m*. 


This can readily accommodate a 10kg/m* change of the density during a run. 


Each density reading of the brine in the upper compartment was acquired by closing 
the vent opening after it was opened for a certain amount of time. The overflow salt 
water collected from the notch was measured and poured back into the upper compart- 
ment. The brine solution was mechanically stirred and then the density was carefully 
measured by the submergence of the hydrometer. The stratification in the upper com- 
partment was corrected by stirring while the lower region was made uniform by the 
magnetic stirring mixer. It was assumed that they were both fully mixed, so the den- 
sity in the lower compartment could be determined from the density in the upper one 


by employing mass conservation. 
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2.2.3 Pressure Measurement 


Since the pressure difference AP between the two compartments is very small and since 
salt water is incompatible with most pressure measurement instruments on the market, 
a considerable effort was made to obtain a suitable pressure measurement device. An 
Omega brand(Model PX154-001DI) low differential pressure transducer was employed 
in the experiment for the sake of its reasonable cost and high accuracy. The resolution 
of the transducer is as low as 0.1% of the full scale which is one volt for zero pressure 
difference and five volts for a one inch water high pressure difference. Since it is only 
compatible with air, it could not be directly connect to salt water, this could damage 
the diaphragm of the transducer. As shown in Fig.2.5 and Fig.2.6, the high pressure 
end of the transducer was connected to tap3, which is the one in the lower compartment 
and is the one nearest to the partition plate, that is, the closest to the opening where 
the flow exchange occurs. Instead of directly connecting the water to the transducer, 


an air column was left at the top of the water column. 


Since in the experiment, the brine level in the upper compartment is constant due 
to the overflow at the notch, the pressure in the upper compartment is constant except 
for the hydrostatic difference. This was verified later by the manometer readings in the 
upper compartment, which indicate that only hydrostatic pressure variation exists in 
the upper region with liquid columns in the manometers being the same as the brine 
level in the upper compartment. So for simplicity, only the high pressure end of the 
transducer was connected to tap3 and the lower pressure end of the transducer was 
left open to the atmosphere. Thus, the pressure value obtained from the transducer 
was the lower compartment’s relative pressure to the atmosphere at location tap3. The 
pressure in the other locations of the lower compartment can be deduced from this 


value by employing hydrostatics. 


Because of the air column and position of the transducer being higher than that of 


tap3, the reading obtained directly from the transducer is not the real pressure value 
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Figure 2.5: Sketch of the pressure measurement arrangement 
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Figure 2.6: Photographs of the pressure measurement arrangement 
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at tap3(Fig.2.3), actually it is Props — Psor9H1 — PamgH2. Since pam is much less than 
Psort, the last term can be neglect. Anyway, the transducer can measure the pressure 
change due to the depression and expansion of the air column, which is caused by the 
changing of the water column, once it was calibrated. The temperature was essentially 
constant and the thermal expansion effect of air column is negligible. Therefore, the 
exact water column change can be obtained and used to figure out the real pressure 
value by using the static pressure formula. The manometer readings from tap] and 


tap2 helped us to check if the transducer reading was within the right range. 


To calibrate the transducer, the lower compartment manometer readings were set at 
the same value as the upper one by moving the storage tank. This implied that all the 
five manometers readings were aligned with the brine level in the upper compartment. 
At this time, both ends of the transducer were left open to the atmosphere. The output 
was one volt representing zero pressure difference. Next, the tube of the manometer3 
was connected to the high pressure end of the transducer by a quick connection fitting 
while the other end of transducer was still left open to the atmosphere. Because of the 
sudden sealing, there was either vacuum or extra pressure inside tube connection which 
changed the output of the transducer away from one volt. The reading was brought back 
to one volt by turning the zero adjustment button of transducer with a screwdriver. 
Then the storage tank was raised to set the lower compartment manometer readings 
0.0254m (1 inch) higher than the upper one. This can be accomplished by observing 
the other two manometer readings in the lower compartment rather than manometer3, 
because of the inaccuracy caused by the air column. The full span adjustment button of 
the transducer was adjusted with the screwdriver to bring the reading back to five volts, 
if it is not, which represents a one inch water pressure difference (248.92 N/m?). This 
procedure was repeated several times to make sure that the reading of the transducer 
is one volt with zero pressure difference and five volts with one inch water pressure 


difference while the storage tank was moved up and down. 


While doing the experiment, the lower compartment manometer level was brought 
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to the desired value according to the transducer reading. As the stopper was removed, 
the water columns of the lower compartment manometers dropped because of the pres- 
sure loss in the connection between the storage tank and the lower compartment. This 
variation was measured precisely from the chartrecorder and the multimeter which were 
connected to the transducer. Thus the pressure difference over the duration of a run 
was obtained accordingly. Results from experiments performed for the same conditions 


indicated that the repeatability of the pressure measurement was within +1.0 percent. 


2.2.4 Fluid Overflow Measurement 


The volumetric overflow measurement is relatively easier. This was done by using a 
beaker, which was placed underneath the overflow collecting box, to collect the brine 
that overflows the notch in the upper compartment during each experimental run. After 
the running time was over and the vent opening was closed, usually there was still some 
continuous brine overflows because of the surface tension at the notch. The duct of the 
overflow collecting box was closed with a rubber stopper when it was confirmed that 
there was no more brine overflow. The total overflow rate was obtained by measuring 
the amount of brine in the beaker with a measuring cup. This overflow salt water was 
then poured back into the upper compartment before stirring the brine. Once the den- 
sity measurement was carried out, the extra amount of brine was drained out through 


the duct at the bottom of the overflow collecting box by removing the stopper. 
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Chapter 3 


Results and Discussion 


3.1 Preliminary Experiments 


Before discussing the results, the definition of the governing parameters is given below 
with reference to Fig.3.1, for combined natural and forced convection flow through a 
horizontal vent. Qy is the volume flow rate of the fluid from the lower compartment 
to the upper compartment, Qp is the volume flow rate of the brine from the upper 
compartment to the lower compartment, Q; is the volume flow rate of the fresh water 
from the storage tank, whose vertical positioning is used for the pressurization of the 
lower compartment, and Qo is the volume overflow rate of the brine at the overflow 


notch in the upper compartment. Therefore, 


Q0=Q1 (3.1) 


Qn is the net upward volume flow rate across the vent opening and, 


Qn = Qu -Qn0=Q0=Q1 (3.2) 


Py is the static pressure at the elevation of the top surface(upper compartment side) of 
the support plate, Pr is the static pressure at the elevation of the bottom surface(lower 
compartment side) of the support plate, both pressure were measured after the vent 
was opened, and AP = Py — Py. Also, py is the density of the fluid in the upper 
compartment and py is the density of the fluid in the lower compartment. The fluids 
in the two regions are well mixed before opening the vent as well as after the vent is 


closed following an experiment. 
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Figure 3.1: Combined natural and forced convection flow through a horizontal vent. 
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3.1.1 Experiments for AP = 0. 


Some preliminary experiments have been carried out for the circumstance of zero pres- 
sure difference across the vent. The upper compartment was filled with brine and the 
lower compartment with pure water. The control valve which connects the storage tank 
and the lower compartment together being kept closed. Thus, no external pressure is 
applied to the lower compartment. The length to diameter ratio L/D of the vent was 
varied, as was the initial brine density in the upper compartment. The exchange vol- 
ume flow rate Q, which is the volume flow rate at zero AP, was thus measured for a 
range of L/D and Ap. Obviously, there is a hydrostatic pressure difference across the 
opening support plate due to the plate thickness and the density difference across the 
plate. Here, AP = 0 means that the external imposed pressure difference is equal to 


zero. Since there is no overflow Qo, 


Qo =91=Qn=0 (3.3) 
Thus, 
Q=Quv=QpD (3.4) 


according to Eqn.(3.2). 


Photographs of the flow exchange phenomenon obtained from the shadowgraph are 
shown in Figures 3.2(a) and 3.2(b) for two different vent openings. A vigorous flow 
exchange was observed to occur in a cylindrical region with its axis coinciding with the 
axis of the vent opening. The brine enters the lower compartment, seen as a plume flow 
region in the lower half portion of the picture, due to the buoyancy force. Fresh water 
enters the upper compartment, seen as a plum flow region in the upper half portion of 
the picture, due to the mass conservation. The flow exchange looks similar to a jet in- 
jection circumstance, but the speed is much lower, since it is due to natural convection 


arising from the density difference between the two compartments. 


The volume exchange flow rate (in m*s~) from the upper compartment to the lower 


compartment, or vice versa, in terms of the rate of dilution of the brine can be shown 


0 = dV WHA szuourtiadxe ayy 10y syderBZojoyd ydeiZmopeys :7'¢ emai 
OT = G/F ‘(%) ¥920'0 = a OT = G/T ‘(™) s¢¥0'0 = a 
‘(gut/64) §L00T = O70 ‘(,us/5y) G'1601 = OHO (q) ‘(gt4/5¥) 96001 = O'Fd ‘(,us/64) 7 6H0T = OH (2) 
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to be given by(see appendix A): 


or —Valden/dt) (3.5) 
(px — pL.o) — y= (PH,0 — pH) 

where Vz; and Vz are the volume of the heavier liquid (brine) in the upper compartment 

and the volume of the less dense liquid (water) in the lower compartment, respectively, 

px is the density of the brine solution at time ft, and pzo and py are the densities of 


the contents of the water- and brine-filled compartments at zero time. Thus, to obtain 


Q, pH was plotted against time t and derivative obtained. 


Two typical plots of the brine density variation with time, with density measured 
by the hydrometers, as mentioned earlier, are shown in Fig.3.3 and 3.4, respectively. 
In each figure, two experiments were carried out with the same L/D ratio but different 
initial brine density pz. The total duration of the experiment (or the total number 
of individual density measurements) was selected to obtain the brine density variation 
as a linear function of time. Best-fit curves were obtained for these linear segments 
of the density-time plots, and the volume rate of exchange flow were thus calculated 
from Equ. (3.5). Note that each density measurement represents opening the vent for 
a specified time interval, mixing the fluid region after closing the vent and measuring 


the resulting density. 


As can be seen from the figures, for the same L/D ratio, the increase of the initial 
brine density px,o in the upper compartment gives rise to an increase in the absolute 
value of the slope dpy/dt. This is reasonable, because higher the py, the stronger is 
the buoyancy force. Thus, the dilution of the brine in the upper compartment, which 
is determined by the exchange flow rate Q and is reflected in the rate of decrease in 
px, is also faster at larger initial density difference. According to Eqn.(3.5), the ex- 


change volume flow rate Q is larger if the rate of decrease in py concentration is greater. 


In Figures 3.5 and Fig.3.6, the volume exchange flow rate Q is plotted in the form 
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Figure 3.3: Brine density versus time for AP = 0 and L/D = 1.0. 
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Figure 3.4: Brine density versus time for AP = 0 and L/D = 0.072. 
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/ 


of a dimensionless Froude number, Fr = Q/(gApD*/p)1/?, as a function of the length- 
to-diameter ratio of the opening L/D for two different sets of experiments. For the 
first set of experiments, the opening diameter ranged from 0.032m to 0.048m and the 
density difference ratio from 0.03 to 0.09. While for the second set of experiments, the 
opening diameter chosen was 0.0127m and the density difference ratio was from 0.05 
to 0.08. Actually, with increasing time t, the brine is diluted because of the flow of 
the water from the lower compartment, and the density in the lower compartment gets 
higher because of the flow of the brine from the upper compartment. Thus, the density 
difference Ap = py — pr across the vent decreases with increasing time ft. So, in the 
calculation of the Froude number Ap is taken as the averaged value over the duration 
of the entire experiment. Also plotted in Fig.3.5 and Fig.3.6 are the flow exchange data 
obtained by Epstein[19] with brine/water as the working fluid. The range of the opening 
diameter and density difference ratio of his experiment was 0.025m to 0.048m and 0.03 
to 0.17, respectively. The experimental conditions and the measured exchange flow 


rates Q are given in Table 3.1. 


For the first set of experiments, the range of the opening diameter D and the average 
density difference ratio Ap/p are within the corresponding range of earlier studies. As 
shown in Fig.3.5., a comparison with the earlier work on this circumstance indicates 
fairly good agreement. The dependence of Q on L/D follows the trend predicted by 
analysis in the earlier study, lending support to these measurements. The experimental 
results were found to be repeatable to within + 1.0 percent, indicating accurate mea- 


surement of the flow rate Q by the chosen procedure. 


For the second set of experiments, the opening diameter employed was much smaller 
than those used in the earlier work. Interestingly enough, in Fig.3.6, it is seen that much 
higher Froude numbers arise for smaller diameter vents than those obtained in the ear- 
lier work with larger diameters. However, the trend of the dependence of Q on L/D 
is similar to the one obtained in the earlier study. This difference for small diameters 


may be due to the way the dimensionless exchange flow rate, Fr = Q/(gApD*/p)*/?, 
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Figure 3.5: Experimental results for countercurrent exchange flow through a horizontal 
vent for AP = 0 (data set I). 
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Figure 3.6: Experimental results for countercurrent exchange flow through a horizontal 
vent for AP = 0 and very small vent diameter, D = 0.0127 m (data set I). 
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Table 3.1: Measured Buoyancy-Driven Exchange Flow Through a Horizontal Vent 
(AP = 0) 


3.43x107 

3.43x1078 
3.18x10.- 
3.18x10- 


0.0445 
0.0445 
0.0889 
0.1333 
0.0159 
0.0159 
0.0254 
0.0254 
0.0350 
0.0350 
0.0508 
0.0508 
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is calculated. Though the numerator Q decreases because of the smaller diameter, the 
D°/? in the denominator decreases even faster leading to higher Froude numbers. This 
implies that the nondimensionalization for very small diameter vents may have to be 


done differently. Also surface tension effects may become significant at these diameters. 


Since the current work is largely directed at the flow exchange with nonzero ex- 
ternal pressure difference AP across the vent, only limited data with AP = 0 were 
obtained. The range of the L/D ratio and the density difference Ap investigated here 
may not appear to be wide enough to establish the generality of Figure 3.6. However, 
results from experiments performed for the same conditions indicated that Fr was re- 


producible to within +1.0 percent leading support to the accuracy of the data obtained. 


3.1.2 Experiments for Unstratified Situation, Ap = 0 


Some preliminary experiments have been carried out for the unstratified circumstance, 
Ap = 0. The upper and lower compartment were both filled with pure water, and the 
pressure was applied to the lower compartment by raising the storage tank. There is 
no buoyancy exchange flow but an upward forced flow Qy arises due to the pressure. 


Also, Q p is equal to zero in this case. Therefore, 


On OT Qo = Qu (3.6) 


The volume overflow rate Qo was measured for different values of the pressure differ- 
ence AP across the vent opening. Here AP is the pressure head maintained by the 
storage tank. The dependence of Qo on AP is shown in Fig.3.7(a). It is obvious that 
an increase in AP gives rise to an increase in Qo, as expected. This variation of Qo 
with AP can be used to develop a correlating equation which will allow an exhibition 
of velocity across the vent at large pressure difference. This is the standard vent flow 
model which was discussed earlier and which is based on Bernoulli’s equation. For 
forced flow through the vent opening due to a pressure difference, Bernoulli’s equation 


gives an estimate for the velocity across the vent as: [2(/AP — gpwaterAH)/pwater}!/?. 
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Figure 3.7(a): Experimental results for the unstratified situation, Ap = 0: Volume flow 
rate Qo versus AP. 
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Figure 3.7(b): Experimental results for the unstratified situation, Ap = 0: Veloci- 


ties obtained by two different methods versus AP. 
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Where Pwater is the density of the pure water, AH the thickness of the opening support 
plate. Also, from the volume overflow rate Qo obtained, average velocity across the 
vent opening is calculated as: Qo/Ay, where Ay is the vent area. Velocities calculated 
from the Bernoulli’s equation and from the volume overflow rate Qo are plotted in 
Fig.3.7(b). As shown in the figure, the former one are substantially higher than the 
later one, since Bernoulli’s equation is an ideal equation, it doesn’t count the pressure 
loss and any other loss caused by the disturbance of the flow. In fact, a flow coefficient 
Cp should be taken into consideration, as described by the standard vent flow model 
[Eqns.(1.3)~(1.5)]. 


3.2 Purging Flow Experiments 


Conditions for unidirectional flow between the spaces connected by a horizontal vent 
opening had been determined by purging flow experiments. The upper and lower com- 
partment contain brine solution and fresh water, respectively. Each experimental run 
involved a unidirectional upward flow of fresh water, from the lower compartment, at 
measured pressure difference AP across the opening and volume overflow rate Qo at 
the notch of the upper compartment. Here AP = Py, — Py is the static pressure differ- 
ence across the support plate for the vent opening. Therefore it includes the hydrostatic 


pressure difference and the pressure difference due to the imposed external force. 


In each experiment, fresh water flows into the lower compartment through the con- 
trol valve at the bottom of the lower compartment due to the pressure head maintained 
by the storage tank(see Fig.2.1). An equivalent volumetric flow of the brine overflows 
from the upper compartment across the overflow notch. At small AP, a limited amount 
of brine solution Qp was observed to flow downward through the opening against the 
upward fresh water flow. In this case, the AP was not strong enough to completely 
overcome the buoyancy force due to the density difference Ap across the opening. At 
higher AP, the brine solution was prevented fon flowing across the opening. This 


situation was evident by the absence of any sign of brine solution below the opening, 
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either as a plume or in the form of individual packets(drops) of fluid. The shadowgraph 
technique was applied here for the observation of these transient pressure. The flooding 
or purging condition was determined as a AP value that represented a transition from 


countercurrent, or bidirectional flow to a unidirectional flow. 


For openings of small L/D, the purging condition was detected by adjusting the 
pressure head of the storage tank so that the downward flow rate of the brine solution 
was reduced to a very slow rate of release of brine fluid packets into the lower compart- 
ment, ultimately reducing to zero downward flow. With the tube opening (high L/D) 
it was possible to adjust the AP so that the brine - water mixing zone filled the tube 


without any net downflow of brine solution. 


The pressure difference AP that produced this condition of zero downward buoyancy 
driven flow was defined and recorded as the purging (or flooding) pressure APrroop, 
and the unidirectional upflow rate of the fresh water Qy was therefore the flooding flow 


rate. Thus, since Qp = 0, we have: 


Qrioop = Qu = Q1= Qo (3.7) 


In Fig.3.8, shadowgraph pictures taken from the purging flow experiments show 
the transition from a bi-directional flow to a uni-directional flow. As AP approaches 
the value of APrroop, but with the pressure difference lower than the purging value, 
tiny down flow of the brine solution was observed as a plume in the lower compart- 
ment(Fig.3.8a). As AP was maintained at APrroop or little higher, a uni-directional 


upflow of the fresh water was found to arise with no observable downflow. (Fig.3.8b). 


The measured flooding pressure APrroop, flooding flow rate Qrzoop and the ex- 
perimental variables are presented in Table 3.2. The dependence of APrrogop on the 


buoyancy force gApD for two different vent opening diameters and various L/D ratios 
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(a) D = 0.0127, L/D = 2.75, px.o = 1062.5(kg/m), AP = 45.92(N/m?) 


(b) = 0.0127, L/D —=2.75) PHO = 1062.5(kg/m?), AR=APrroop = 53.70(N/m7?) 


Figure 3.8: Shadowgraph pictures of purging flow experiments: Uni-directional up flow. 
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Table 3.2: Measured Uni-directional Purging Flow Through a Vent Opening, while 
py = 1000.0kg/m*. 


Qrroop(m's"*) x 10-° 
1.27x10- ; | 
1.27x10-5 

127105 

1.27x10-3 

1.27x10-3 


0.0159 
0.0159 
0.0159 
0.0159 
0.0159 


0.0762 
0.0762 
0.0762 
0.0762 
0.0762 


0.0508 
0.0508 
0.0508 
0.0508 
0.0508 
0.0508 


0.0350 
0.0350 
0.0350 
0.0350 
0.0350 


0.0127 
0.0127 
0.0127 
0.0127 
0.0127 


3.18x10-? 
3.18x1073 
3.18x10-3 
3.18x10-3 
3.18x1073 
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are shown in Fig.3.9(a) and Fig.3.9(b), respectively. As seen in Fig.3.9, the purging 
pressure APryroop increases with an increase in gApD, which represents the buoyancy 
force due to a density difference Ap across the opening whose diameter D is taken as 
the characteristic dimension. This is reasonable, since higher the buoyancy force across 
the opening, higher the purging pressure needed to overcome this force to give rise toa 
uni-directional upflow. It is found that the dependence of APrroop on gApD follows 


a square root curve behavior. 


The opening length to diameter ratio L/D does have an effect on the purging pres- 
sure APrroop. As seen from Fig.3.9(a), an L/D value of around 1.0 is the turning 
point. As L/D decreases below this value to L/D = 0.25, the purging pressure dif- 
ference APrroop also decreases. As L/D increases from this value to L/D = 6.0, an 
opposite behavior is noted, with the purging pressure APrzoop decreasing with in- 
creasing L/D. There is an intersection of the two pressure curves for L/D = 1.0 and 
L/D = 2.754. At smaller density difference Ap, APrroop increases with increasing 
L/D. At larger Ap, APrzoop is found to decrease with increasing L/D. At extremely 
high L/D ratio, which is a tube opening, the purging pressure APrroop needed to 
generate the uni-directional upflow is much smaller than that at moderate L/D values. 
As L/D goes to the other extreme, the orifice, which has a very low L/D ratio, the 
purging pressure APrroop needed to generate the uni-directional upflow drops down 
to a lower value again. Clearly, L/D is an important parameter and affects the flow 
patterns and regimes that arise, as discussed by Epstein[19]. A more detailed effort on 


the flow field is therefore needed to fully understand these effects. 


The opening diameter D also has effect on the purging pressure difference APrroop. 
As tabulated in Table 3.3 and depicted in Fig. 3.10, when D increases from D = 0.0127 
m to D = 0.0159 m, APrroop becomes higher for the same L/D ratio and density 
difference Ap range. So, for a larger opening diameter, the purging pressure difference 
APrroop must be maintained at a higher value to prevent any downflow of brine so- 


lution from entering the lower compartment. A higher opening diameter D leads to 


A Prroop, N/m 


A Prtoop, N/m 
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Figure 3.9: Purging flow experiments: The dependence of APrroop on gApD 
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Figure 3.10: Purging flow experiments: The effect of the opening diameter D on purging 
pressure. 
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a stronger buoyancy force gApD. Thus, the purging pressure needed is also higher. 
Therefore, the purging pressure depends on the diameter D through the buoyancy force 
and an L/D which represents the geometry of the vent and which affects the flow pat- 


tern. 


Considering now the flooding flow rate Qrzoop, it is found that the dimension- 
less Froude number, Frppoop = ae. does depend on the density difference 
Ap, along with its dependence on the L/D ratio. Actually, as tabulated in Table 3.2, 
QFrroop is relatively scattered for a given L/D ratio but different Ap, since Orzoop 
is determined by the two other parameters: APrroop and the density difference Ap 
across the opening. When APrzoop is maintained at a very high value, the density 
difference Ap, which reflects the buoyancy force to be overcome by the purging pres- 
sure APrroop, is also very high. Thus, the resulting Qrzoop is found to be lower. 
However, Frpioop is averaged over different Ap for a given L/D ratio and plotted in 
Fig.3.11, from which the dependence of Frrioop on L/D is shown to be similar to the 
one obtained by Epstein and Kenton[22](see Fig.1.22). The magnitude of the Froude 
number obtained here is relatively higher. The reason is the same as the one in section 
3.1.1, because the opening diameters employed are much smaller than those in the early 
work. In fact, Epstein’s conclusion of the independence of F'rppoop on Ap was based 
on his former experiment with purely buoyancy force, only the density difference Ap 
determine the exchange flow rate for a given opening size, and the dimensionless flow 
rate, Fr = Q/(gApD*/p)*/?, eliminates the effectiveness of the Ap. But for the forced 
upward flow with buoyancy force, as mentioned before, pressure difference AP also 
plays a role in Qrzoop, and it is not reflected in dimensionless Froude number. Thus, 


more scatter data was obtained. 


The dependence of the dimensionless Froude number F'rp,oop on the dimensionless 
parameter AP/gApD, for different Ap across the vent opening is shown in Fig.3.12. 
The data fall well within a band, indicating that the increase AP/gApD gives rise to 


an increase in Frpyoop. Here, AP/gApD is a buoyancy parameter which reflects the 
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Figure 3.11: Purging flow experiments: Dimensionless purging flow rate versus opening 
aspect ratio L/D. 


62 


nN 
. D=0.0127m , 

M= Ap = 147.4 kg/m, 

O= Ap = 109.4 kg/m, ° 

@= Ap = 90.08 kg/m, nd 
° O= Ap = 63.90 kg/m 
r, j 

4 


Fretoop 


1.8 2.1 2.4 2.7 3.0 


0.0 0.3 0.6 0.9 ice 1.5 
A Prroop/g4pD 


Figure 3.12: Purging flow experiment: The dependence of the purging flow rate on 
dimensionless parameter AP/gApD. 


63 


tatio of the external force to the buoyancy force. Higher the gApD, the higher the 
AP needed to give rise to a given purging flow rate. Thus, this nondimensionalization 


eliminates the characteristics of different purging flows. 


3.3. Combined Flows 


The main intention of this portion of the study is to determine the effect of an imposed 
pressure difference on the mass flow rate across a vent opening in a horizontal partition. 
This flow exchange through a horizontal vent involves a combination of a buoyancy 
driven exchange flow and a forced flow. Buoyancy force due to the density difference 
across the vent opening gives rise to a bi-directional exchange flow through the vent. 
The imposed pressure difference gives rise to a forced upward flow through the vent. 
Thus, a mixed natural and forced convection flow arises due to the combination of these 
two forces. The apparatus used for the combined or mixed convection experiments is 
the same as the one for the purging flow experiments. Combined convection conditions 
were obtained simply by adjusting the value of AP, either raising or lowering the storage 
tank, to values below the corresponding purging pressure difference APrroop. The 
rate of dilution of the brine in the upper compartment was obtained by measuring the 
density and the volume flow rate of the overflow brine Qo at the notch of the upper 
compartment. The buoyancy driven downflow of the brine solution Q p and the forced 


upflow of the water Qy can be shown to be given by (see appendix B): 


—Vudpy/dt 
pay ee hed SSF et 3.8 
peeves rE (38) 
—Vidpy/dt 
apy Sanaa cay (3.9) 
PH — PL 
and 
V; At 
PL = pL — —(ox — pHo) — iy: — Pwater) (3.10) 
Vz Vz 


Where, Vy and V1, px and pz, px. and pro are the same as those defined for Eqn.(3.5), 
Q1 is the volume flow rate from the storage tank to the lower compartment, giving 


Qr1 = Qo, At is the time interval from the starting point to the current point, and 
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Pwater is the density of the fresh water entering from the storage tank to the lower 
compartment. The density of the water in the storage tank is that of pure water at 


room temperature, that is about 1000.0 kg/m°. 


The shadowgraph photographs of combined natural and forced convection flow are 
shown in Figs.3.13 to 3.14. Some are with the down flow Qp being smaller than 
upflow Qy and some are with down flow Qp being almost the same as the upflow Qy, 
depending on the initial density difference and the pressure difference maintained. It 
was observed that the exchange phenomenon is similar to that with AP = 0, except 
that amount of flow in the two directions need not to be equal now. For the experi- 
mental conditions chosen, most of the cases are with Qy greater than Qp. Thus, the 


effect of the imposed pressure is greater than that due to the buoyancy force. 


The measured external pressure difference AP, the net volume overflow rate of brine 
at the notch of the upper compartment and the experimental variables are presented 
in Table 3.3. Also given in table 3.3 are the buoyancy driven component Qp and the 


forced flow component Qy of the flow rate determined from the experiments. 


3.3.1 Density Variation Versus Time 


Four typical plots of the brine density variation with time, as obtained from the hydrom- 
eter measurements, are shown in Figs.3.15(a) Tega (d). In each figure, two different 
experimental results, which were carried out for a given L/D ratio and essentially the 
same initial brine density py but different pressure difference AP, are plotted to de- 
pict the effect of the externally imposed peeeere on the combined convection transport. 
A linear best fit analysis was carried out for the linear segments of the density-time 
plots. It can be seen that an increasing AP value gives rise to an increasing absolute 
value of the slope dpy/dt. This, in turn, indicates an increase in the forced upward flow 
component Qy, from the lower compartment to the upper compartment, according to 
Eqn.(3.8). This is also confirmed by the calculations predicted in the following sections. 


The buoyancy driven downward flow component Qp, from the upper compartment to 
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(a) D = 0.0127 (m), L/D = 2.75, py = 1050.5(kg/m?), pr = 1001.0(kg/m*), 


AP = 61.49 (N/m?) 


(b) D = 0.0127 (m), L/D = 2.75, px = 1050.5(kg/m?), pr = 1000.0(kg/m?), 


AP = 52.77 (N/m?) 


Figure 3.13: Combined flow under buoyancy and the imposed pressure: Shadowgraph 
photographs 
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(a) D = 0.0127 (m), L/D = 2.75, px = 1090.6(kg/m°), pr = 1000.0(kg/m’), 


AP = 36.51 (N/m?) 


(b) D = 0.0127 (m), L/D = 2.75, px = 1089.0(kg/m*), pr = 1000.5(kg/m?), 


AP = 39.95 (N/m?) 


Figure 3.14: Combined flow under buoyancy and the imposed pressure: Shadowgraph 
photographs 


Table 3.3: Measured Combined Convective Flow Through a Vent Opening 
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0.0763 


0.0508 
0.0508 
0.0508 
0.0508 
0.0508 
0.0508 


0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 
0.0350 


0.0127 
0.0127 
0.0127 
0.0127 
0.0127 
0.0159 
0.0159 
0.0254 
0.0254 
0.0254 
0.0254 
0.0254 
0.0254 


3.18x10-5 
3.18x10-3 
3.18x10-$ 
3.18x10-% 
3.18x10-5 
3.18x10-3 
3.18x1075 
3.18x10-3 
3.18x10-3 


67 


68 


‘oulty 


SNSIOA UOTeIIeA AIsUa(] :MOY UOTZI9AU0D padIOJ pue [eMyeUu pauIquUIOD :cT'¢ emM3Iy 


“upus ‘auagy 
ow oor ee oe ec oe oe oe oc om oy ee 
3 
3 
g 
3 
sf 
33 
o 
= @ 
38 
se 
(p) =o 
33, 
3 
: 3 
(= /Y) ¢-C1* 916 F = P/*4P 1, /N) Mm aT 3 
(y= /P4) gO LIN © He/"4p '(,e/m) te eta 
v= a/1 (=) wI00-0 - 
F 
“UyU “aug 
ow oe oe oe ee oe os er oc oe oe ee 
he 
2. 3 
: 
| 3 
eo 
| : 
os wo 
et 
oe BS 
(4) af 
oe 
og 
AS 
| bd » 


ose 


Ooo 


(= /94) y.01- 1969 WN) 26 tt~ag 


diteo=a 


Oarr 


“UyUn “aut 


ost oes om co ee we oo oe ec Lal ee 


3 
z 
e 
e 
te 
3 
$9 
°s 
3¢ 
i 
35° 
8 
cs 
(42/79) «01% 9069 = 19/847 Lin) mea = 8 
(y= /PQ) «1 OLR: = 19/847 “= /n) tL 8 
e-a/7 l=) wise=g = 
s 
2 
“Uyus ‘oun 
oe oo on oa oo os ee oe os oe 


cro! Tere! 


Varor 


o 
EH 
Fe] 
3 @ 
28 
e353 
_ 
3% 
or 
°s 
.3, 
J 3 
(s,=/PQ) g-OF° OLE = 1P/94P (= /m) tk OMA e 
(9, /P4) y-01 «C82 & = 1P/*4P “|, /N) OC Ota 
e- 0/1 (=) ieee = 
2 
: 
= 


(>) 


69 


the lower compartment, is determined by both Qy and Qo and is difficult to predict 


in terms of the governing variables. However, the basics trends are discussed below. 


3.3.2 Volume Flow Rate Results 


Typical plots for the volume over flow rate Qo which is equal to the net volume flow rate 
Qn, versus Ap, at fixed values of AP, opening diameter D and L/D ratio, are plotted 
in Figs.3.16 through 3.19. Also plotted in these figures are the buoyancy driven down- 
ward flow Qp and the forced upward flow Qy. The flow rates Qo and Qy were found 
to decrease with an increasing density difference Ap, while Qp was found to increase. 

The reason for this behavior is the same as that for the purging flow experiments. 
Higher the value of Ap, stronger is the downward buoyancy force, leading to a larger 
downward flow Q p, and a weaker forced upward flow Qy for a given pressure difference 
AP, since larger pressure is needed to overcome the buoyancy force. This also results in 
a smaller net overflow rate Qo according to Eqn.(3.2). Therefore, the density difference 
Ap across the vent opening has a negative effect on the volume overflow rate Qo and 


forced upward flow Qy, but a positive effect on the buoyancy driven downward flow Q p. 


Typical plots for the net volume flow rate Qo or Qn versus AP, at given values of 
Ap, opening diameter D and L/D ratio, are plotted in Figs.3.20 to Fig.3.23. Also 
plotted in Fig.3.22 and Fig.3.23 is the forced upward flow component Qy. Qo was 
observed earlier to increase with an increasing AP in the experiment. Qy was found 
to increase with an increasing AP according to the calculation, while Q p was found to 
fluctuate as AP was increased and did not give out any straightforward dependence on 
AP. Thus, Qp is not shown in the figures. Consequently, the imposed pressure differ- 
ence AP across the vent opening has a positive effect on the volume overflow rate Qo 
and forced upflow rate Qy. This is because the imposed positive pressure difference AP 
was designed to oppose the downward buoyancy force which is caused by the density 
difference Ap across the vent opening. It is also noted that, for given opening diameter 


D and L/D ratio, AP needs to be maintained at higher values for experiments with 
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Figure 3.16: Combined forced and natural convection flow: Volume flow rate Qo versus 
density difference Ap at AP = 42.75N/m?, D = 0.0127 panty Vey PS em (ay waa 
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Figure 3.17: Combined forced and natural convection flow: Volume flow rate Qo versus 
density difference Ap at AP = 80.20N/m?, D = 0.0254 m, L/D = 1.0 
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Figure 3.18: Combined forced and natural convection flow: Volume flow rate Qo versus 
density difference Ap at AP = 57.18N/m?, D = 0.0127 m, L/D = 1.0 
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Figure 3.19: Combined forced and natural convection flow: Volume flow rate Qo versus 
density difference Ap at AP = 71.90N/m?, D = 0.0159 m, L/D = 1.0 
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Figure 3.20: Combined forced and natural convection flow: Volume flow rate Qo versus 
pressure difference AP at fixed Ap, D = 0.0127, L/D = 0.25. 
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Figure 3.21: Combined forced and natural convection flow: Volume flow rate Qo versus 
pressure difference AP at fixed Ap, D = 0.0127, L/D = 2.754. 
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Figure 3.22: Combined forced and natural convection flow: Volume flow rates Qo and 
Qu versus pressure difference AP at Ap = 84.42kg/m?, D = 0.0254 and L/D = 1.00. 
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Figure 3.23: Combined forced and natural convection flow: Volume flow rates Qo and 
Qu versus pressure difference AP at Ap = 106.58kg/m°, D = 0.0127 and L/D = 4.00. 
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higher Ap, in order to obtain the same amount of net volume flow rate Qj compared 
to the experiments with lower Ap(Fig.3.20 and Fig.3.21). This is obviously because a 
larger external pressure is needed to balance the downward buoyancy force for a higher 


density difference. 


In Fig.3.24 and Fig.3.25, the volume flow rate Qo is plotted in the form of the di- 
mensionless Froude number, Fro, as a function of the dimensionless mixed convection 
parameter AP/gApD. It is seen that Fro increases with an increase in AP/gApD. 


This is again an expected behavior as discussed above. 


3.4 Correlations for Combined Flow Through a Horizontal Opening 


Correlation equations have also been derived for the combined flows, based on the 
experiments of uni-directional upward flows and bi-directional combined natural and 
force convection flows. These equations were developed to investigate the dependence 
of the net volume flow rate of Qn , which is equal to volume overflow rate Qo, on 
parameters such as the pressure difference AP across the vent opening, the density 
difference Ap across the opening and the L/D ratio. The equation obtained from the 


data is shown below as: 

Qo = 61x 10m (AP) a (Ap) moa Doe: (3.11) 
The parameter ranges covered by this equation are: 

0.25 < L/D < 6.08, and0.043 < Ap/p < 0.13 (3.12) 


The correlation coefficient for the equation obtained is 20 percent. The above corre- 
lating equation indicates that Qo has a strong dependence on the imposed pressure 
difference AP, since it is proportional to AP to the power of 0.83, it shows that the 
pressure difference has a positive effect on Qo which verifies the conclusion obtained 
in previous section. Qo varies with the density difference Ap across the opening as a 


power of -0.29, the effect of Ap on Qo is little weaker comparing to the one of AP, and 
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Figure 3.24: Combined forced and natural convection flow: Dimensionless volume flow 
rate F’ro versus dimensionless mixed convection parameter AP/gApD. 
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’ Figure 3.25: Combined forced and natural convection flow: Dimensionless volume flow 
tate F’'ro versus dimensionless mixed convection parameter AP/gApD. 
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Ap has an negative effect on Qo since Ap is in the denominator, which also verifies the 
conclusion obtained by previous section. The variation of the opening ratio L/D has 
a small effect on the volume flow rate Qy. The correlation shows that over the L/D 
range chosen, Qy varies with L/D as a power of 0.11. Actually, at smaller L/D values, 
Qo is proportional to L/D, but at larger L/D a negative exponent is obtained. Thus, 


“0.11” is the overall exponent for the whole L/D range. 


A similar correlation equation had been obtained for the dependence of the uni- 
directional upward flow Qy on the pressure difference AP, the density difference Ap 


and the opening ratio L/D. The equation is shown as: 

Qu = 5.01 x 10-8(AP)*>(Ap)-°(L/D)°* (3.13) 
where the parameter ranges covered are: 

0.25 < L/D < 6.08, and0.043 < Ap/p < 0.13 (3.14) 


The correlation coefficient for the equation obtained is 18 percent. The equation shows 
that the effects of AP, Ap and L/D on Qy are analogous to those on Qo, except that 
the density difference Ap has stronger influence on Qy and the pressure difference AP 
has weaker influence on Qy, comparing to those on Qo. For the effect of the density 
difference Ap, since Qy is forced upflow, it is largely generated by imposed pressure 
difference which overcomes the downward buoyancy force due to the density difference 
Ap across the vent. In other word, more downward buoyancy force needs to be bal- 
anced in Qy than in Qo. For the effect of the pressure difference AP, since Qo is the 
difference between the volume upward flow rate Qy and the volume downward flow 
rateQ p, which are both determined by the pressure difference AP. Thus, the pressure 


difference AP has stronger effect on Qo than on Qu. 
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Chapter 4 


Conclusion and Future Work 


4.1 Summary and Conclusions 


A series of experiments has been carried out to investigate the combined natural con- 
vection and forced flow through a single horizontal] vent opening in the ceiling of a 
room containing a fire. The buoyancy driven flow, without an imposed pressure, and 
the forced flow, without buoyancy effects, were also examined. A plexiglass tank with 
a vented partition in the middle was constructed to simulate the room with a fire and 
the ambient environment, by filling the upper and lower compartments with brine and 
pure water, respectively. Thus, the density difference is obtained due to the salt con- 
centration. The lower compartment could also be raised to a pressure higher than the 
ambient by providing an inflow from an external tank of pure water that could be raised 


or lowered to vary the pressure. 


The average net volume flow rate Qo and forced volume upflow rate Quy were ex- 
perimentally obtained for a range of variation in the pressure difference AP across the 
opening, density difference Ap across the opening and the opening aspect ratio L/D. 
For each geometrical configuration (fixed L/D ratio), the pressure difference across the 
opening was varied, while the density difference was maintained at a certain fixed value, 
and vise versa. Flow visualization using the shadowgraph technique was employed to 
study the flow pattern in the vicinity of the vent opening. Data on the mean volume 
flow rates Qy, Qo and Qp were obtained by measuring the brine density in the upper 


compartment resulting from the flow exchange across the opening. 


The experiments demonstrated that for combined natural and forced convection flow 
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and for the range of opening ratio L/D, density difference Ap and pressure difference 
AP considered, the pressure difference has the strongest effect on the volume flow rates 
Qo and Qy. Both of these were found to increase with an increase in the pressure 
difference. The density difference Ap across the opening, however, has a negative effect 
on the volume flow rates Qo and Quy, the effects being weaker compared to the effect of 
the pressure difference, over the ranges considered. The opening ratio L/D was found 
to have only a small effect on the volume flow rates. In the smaller L/D range, an in- 
crease in the L/D ratio gives rise to an increase in the volume flow rates. In the larger 
L/D range, an increase in the L/D ratio gives rise to a decrease in the volume flow 
rates. It was found that, for the range of parameters considered, the average volume 
flow rates could be correlated fairly accuratly using two separate correlations. For the 
different volume flow rates considered, the effect of the same parameter was found to 
be quite different. For example, the effect of the pressure difference on Qo is stronger 
than the one on Qy with the exponent of the pressure difference AP in Eqn.3.13, the 
correlating equation, being higher than the one in Eqn.3.12. Similarly, the effect of 
the density difference on Qo is weaker than the one on Qy, with the exponent of the 
density difference in Eqn.3.11 being lower than that in Eqn.3.14. 


Purging flow experiments were carried out to determine when the downward flow 
due to buoyancy force is completely stopped by the imposed pressure difference. It was 
found that the density difference across the opening and the opening size determine 
the flooding pressure APrroop required to generate a uni-directional upflow. That is, 
higher the density difference or larger the opening diameter, higher is the APrroop 
needed to prevent any buoyancy driven down flow from entering the lower compart- 
ment. It was found that the density difference also affects the dimensionless flooding 
flow rate. The deviation in Frppoop was more than 20 percent as Ap/p varied from 
0.03 to 0.17. This indicates that the dimensionless flow rate used in earlier studies dose 
not correlate the data satisfactory over the entire parameter range. But the trend of 
the variation in the averaged dimensionless flooding rate with the opening L/D ratio 


agreed very well with the previous work. 
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Experimental results on pure natural convection, with no externally imposed pres- 
sure, were compared with a earlier study, indicated that, within the same opening 
diameter range, the results were in very good agreement. For much smaller opening 
diameters, higher dimensionless Froude numbers were obtained at the same opening 
L/D ratio, compared with the previous study. Again, this suggests the need for a 
different nondimensionalization. Experiments with forced uni-directional upflow and 
without buoyancy effects indicate that an increase in the externally imposed pressure 
gives the rise to an increase in the upflow, which is also the net flow across the opening. 
The velocity calculated from the measured pressure difference AP across the vent, us- 
ing Bernoulli’s equation, agrees very well with the velocity obtained from the overflow 
measurement, if a flow coefficient is taken into consideration. This leads support to the 


pressure measurement procedure employed in the current study. 


Flow visualization experiments revealed a fairly symmetric one way or two way 
flow pattern in the vicinity of the vent opening. For a zero pressure difference, a bi- 
directional symmetric flow across the vent was observed, with the amount of the flow in 
upward and downward directions being equal. As the external pressure was increased 
from zero, the amount of the buoyancy-driven down flow was observed to decrease. As 
the external pressure difference reached a certain value, for the parametric experimental 
conditions chosen, the buoyancy-driven down flow was observed to drop to zero, i.e., the 
flow was prevented from entering the lower compartment. Thus, only a uni-directional 


upflow was generated, this being the purging flow. 


4.2 Future Work 


A Dantec brand hot film probe(Model 55A76) and a constant temperature hot wire 
anenometer (Model 55M01) had been procured for velocity measurements near the 
vent opening. A plexiglass positioner was fabricated to move the hot film to the desired 


location. It can be moved over the X-Y plane and the height of the hot film can also be 
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adjusted. Thus, the hot film can be either placed at the top center of the vent opening 
during the experimental run or at the corner of the plexiglass tank during the mixing 
process. Due to time limitations, experiments for the velocity measurements could not 
be carried out. Experiments on the velocity distribution at the top of the vent opening 
(upper compartment side) will provide inputs to an analysis of the transient flow across 


the vent and the pressure distribution in the vicinity of the vent opening. 


Further improvements in the pressure measurements, with a salt-water compatible 
transducer, will eliminate the inaccuracy caused by the air column used in the current 


study. The inconvenience caused by the calibration needed before each experiment can 


also be eliminated. 


Experiments can also be carried out with square horizontal vent openings, with 
different aspect ratio §, defined as the ratio of the opening height to the side of the 
square opening, that is, L/S. The results can be compared with the results with circu- 
lar openings by making the characteristic identification D = 1.15. It is expected that 
these two results would be similar. Also, inclined vents are of interest in many pratical 
circumstances and need further investigation. The effect of multiple openings on the 


flow and transport processes is another area of considerable importance. 
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Appendix A 
Exchange Volume Flow Rate Through Vent Opening 


with Zero External Pressure Difference 


Equation (3.5) was derived on the basis of mass conservation. The mass balances on 


the brine and water compartments give 


d 

Vi-Se = -Qpu + Qo1 (A.1) 
dpy 

ome = -Qpi+Qpx (A.2) 


Solving equation (A.1) for Q yields 


_ —Vudpy/dt 
PH — PL 


Q (A.3) 


Adding equation (A.1) and (A.2) and integrating the result from the initial time, with 


densities pz9 and pzo we obtain 


Vu(ex — pH,0) + Vi(px — pro) = 0 (A.4) 


Equation (3.5) is then derived by solving the above expression for pz and substituting 
the result into the equation (A.3) as: 


—Vu(dex/dt) 


~ (pu — pio) — ¥2(PH0 — pu) oe 
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Appendix B 
Volume Flow Rate Through Vent Opening with 


Non-Zero External Pressure Difference 


Equations (3.8) through (3.10) were obtained in a similar way to that as follows: The 


discussion in Appendix A Mass balances on the brine and water compartments give 


d 
Vee = —Qppx + QueLt — QopH (B.1) 
d 
Vea = —Qupi. + QppeH + Q1Pwater (B.2) 
Also, 
Q1=Qo0=Qu-Qp (B.3) 


Adding Equation (B1) and (B2) and integrating the result from the initial time, with 


densities pzo and pro, we obtain 
Vi(ex — px,0) = —Va( px — pH) — (PH — Pwater)QrAt (B.4) 
Equation (3.10) is then derived by solving the above expression for pz: 
Vi Q,;At 
- sign ae) ¥ pr Aaa, aL Ts. B.5 
PL = PLO ¥,- (pH — pH.) Vi, (PH — Pwater) (B.5) 


Substituting Equation (B3) into the Equation (B1), Equation (3.8) and (3.10) were 


thus obtained: 
—Vydpy/dt 
= B.6 
bere Apiezeo (38) 
—Vyd dt 
One —Vudpu/dt _ Qo (B.7) 


PH — PL 
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